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The Jeans instability is analyzed for dense magnetohydrodynamic plasmas with intrinsic magneti- 
zation, the latter due to collective electron spin effects. Furthermore, effects of electron tunneling as 
well as the Fermi pressure are included. It is found that the intrinsic magnetization of the plasma will 
enhance the Jeans instability, and can significantly modify the structure of the instability spectra. 
Implications and limitations of our results are discussed, as well as possible generalizations. 
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I. INTRODUCTION 

Whenever the internal pressure of a gas or plasma cloud is too weak to balance the self-gravitational force of a mass 
density perturbation, there will be a collapse. The study of such mechanisms was pioneered by Jeans, and it is well 
established that the Jeans instability plays a crucial role in the mechanism responsible for the formation of structures 
in the universe [l[. 

For lengthscales comparable to the dc Broglic wavelength of the charge carrier, tunneling effects become important 
and give rise to dispersion. This is particularly important for low temperature and/or high density plasmas. A 
characteristic condition for this effect to be of significance is when the de Broglie wavelength is comparable or larger 
than the mean separation of particles. The plasma then behaves as a Fermi gas, obeying Fcrmi-Dirac statistics rather 
than Boltzmann-Maxwell used in classical plasmas. These effects can be captured within the magnetohydrodynamic 
(MHD) model through a modified equation of state for the pressure @, 0|, and by including the Bohm potential 
0, HI H, H, 3 describing quantum forces due to non-locality effects. Quantum plasmas was first studied by Pines in 
the 1960's a0]i an d niany studies has appeared since then [13], e.g., kinetic models of the quantum electrodynamical 
properties of nonthermal plasmas and covariant Wi gner function descriptions of relativistic quantum plasmas 
[121 ]. The s tudy of quantum plasmas in recent years have been motivated by e.g. possible applications to nanoscale 
technology (Hf, new developments in microelectronics [HI, the discovery of ultracold plasmas [IE Ql| [TtJ and the 
experimental demonstration (T|| of collective modes in ultra cold plasmas, new laser plasma/solid-matter interaction 
regimes offered by the next generation of high intensity light sources [T3, [13, [U [22j , as well as developments in laser 
fusion 12311 . Quantum effects are also believed to be of importance in the interiors of compact astrophysical objects 
[2~il. I25l.l2q| such as white dwarfs, neutron stars, magnetars and supernovas, where the density can reach ten orders of 
magnitude that of ordinary solids. Furthermore, in some astrophysical environments, such as in the vicinity of pulsars 
[27l I28I] and magnetars [29 ] . strong external magnetic fields may be present. In such dense and/or strong magnetic 
field environments a quantum description of the plasma which incorporates the spin of t he p article is needed. A great 
deal of interest has been directed toward finding such quantum plasma descriptions [13, [H[ HI, HE HE HE HH ■ I n the 
context of the Jeans instability, quantum effects have been studied in Ref . [371 ] . 

The objective of this paper is to investigate the Jeans instability in a magnetized self-gravitating spin-1/2 plasma. 
Starting with the MHD equations, Poisson's equation for the gravitational potential, and an expression for the mag- 
netization of the plasma due to the electron spin, we derive a dispersion relation for plasma modes with arbitrary 
directions of propagation. The overall stability of the system is investigated, and in particular we study the sta- 
bilizing/destabilizing effect of the Bohm potential and the electron spin. Finally, we discuss the implications and 
limitations of our model. 

II. GOVERNING EQUATIONS 

The governing equations for self-gravitating spin- 1/2 ideal MHD plasmas [E HI, HI HI] are the mass density 
conservation law 



d tP + W-(pv) =0, (1) 



the momentum equation 



p($ + v.V)v = -pW-Vp-V( £--m-b) + B- vf— -m) +-^V (?-^!.) , 
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the idealized Ohm's law 

3 t B = Vx(vx B), (3) 

and Poisson's equation for the gravitational potential 

V 2 = AirGp. (4) 

Here p is the mass density, v is the fluid velocity, B is the magnetic held, M is the magnetization due to the electron 
spin, (j> is the gravitational potential determined by Poisson's equation, K is Planck's constant divided by 2tt, m e> i is 
the electron and ion mass respectively, and G is the gravitational constant. The momentum equation, eq. ([2]), has 
been modified to include non-locality effects, such as tunneling, described by the Bohm potential. These effects are 
important for dense and/or low temperature plasma systems and for short wavelengths. 

This system of equations is closed with an equation of state for the pressure and an expression for the magnetization. 
The equation of state for the pressure is written as, VP = c 2 Vp, where c s = (dp/dpo) 1 ^ 2 is the sound speed. The 
ion-acoustic velocity can be written as 



c l = < + ^ ( vl + -v 2 Fe ) , (5) 
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where Vu and Vt e are the ion and electron thermal velocities and v Fe = 7L(37r 2 ne) 1//3 /m e is the electron Fermi velocity. 

1/3 

For temperatures well below the Fermi temperature, the thermal velocity is given by v ts = C(kgT/hn s ), where C 
is a dimensionless constant of order unity, whereas for temperatures much higher than the Fermi velocity we have 
v ts = (IcbT '/rris) 1 / 2 . Here n s is the number density of particle species s, ks is Boltzmann's constant and T is the 
temperature. 

Multifluid equations including the effect of the magnetization due to the electron spin has been derived in Ref. (39j 
and further developed for MHD regimes in Ref. [32j ■ It is then showed that for dynamics on a time scale much slower 
than the spin precession frequency, the magnetization is given by 



M =^tanh(^) 

mi \ k B T J 



B, (6) 



where, B = B/_B, and B = |B|. Here it has been assumed that the spin orientation has reached the thermodynamic 
equilibrium state in response to the magnetic field, which accounts for the t&n}i(pB B / kgT)-ia,ctor , where ps = e?i/2m e 
is the Bohr magneton. Furthermore, on a time scale shorter than the spin relaxation time scale, the individual electron 
spins are conserved, and thus we can take tanh(/xs B / ksT) as constant for an initially homogeneous plasma (40| . 

Assuming perturbations on a homogeneous background (equilibrium values denoted by index 0) we write p = pa+5p, 
= 4>a + 5(f), v = 5v, B = Bqz + (5B, and M = Mqz + <5M and linearize our governing equations. Making a harmonic 
decomposition, we obtain 

ujSp — pq\l ■ 5v = 0, (7) 

- iojp a 6v = -ik (p 8<t> + c 2 s 5p + B ° 6Bz - M Q 6B Z - B SM^\ + ik z B Q (— - Sm] - ^ k k Sp, (8) 
V Mo / V Mo / 4m e m,i 

~ codB = B k x (Sv x z), (9) 
- k 2 5<f> = AnGdp, (10) 



and 



B po 

where <5B_l = (5B X ,SB V ,0), and M = (B /Mo)x/(1 + %)• Here 

^m. tanh ( hbBo 

mid \ Kb -i 



SM = M ( + — Z I , (11) 



1 - ^& tanh ( mZa. 



(12) 
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is the magnetic susceptibility and Ca = (-Bo/MoPo) " is the Alfven velocity. 

Eq. @ already gives us an expression for (5B in terms of 5v and we can solve eq. ([7]) for 5p, and then solve eq. |TT 
for S(j). Using these results in eq. ([8]) we find an equation for 5v. The parallel and perpendicular component of <5v are 
given by 



u> 5v z 



2l2 



K 4 k 



k z k. ■ <5v 



X 



i + x 



C A kzk. ■ 5v± 



and 



lu S\r± = — kj_ 



AnGpo 



k + c 2 k 



1 



2 7,2 



-cik 



k • 5\ - 



1 



(13) 



-C 2 A k 2 z 5w x . (14) 



fc 2 " ' " s " 1 l + x l + x "~ ' 4m e mr7 " 1 + 

Eqs. (|13p and (|14[) can be written in matrix form as D ab Vb = 0, and the dispersion relation is obtained from ||-D ab | = 0, 



1 



1 



-C\ k 2 



X 
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i + x 



c A 



2 k 2 
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Here we have used the notation, 



V 2 



c 2 + 



h z k 



4m P m,- 



(16) 



The first factor of the dispersion relation describes the shear Alfven wave, which is always stable. Focusing on the 
stability properties, this mode is not considered below. The second factor describes the fast and the slow magne- 
tosonic modes modified by a gravitational potential, non-locality effects due to the Bohm potential, and the spin 
magnetization. Our dispersion relation agrees with that found in Refs. [32j and j3o| in the appropriate limits if we 
correct for a sign error in those papers. 

In order to gain some understanding of the stability properties, we first study the special cases of parallel and 
orthogonal propagation to the background magnetic field. Initially we discard the contribution h 2 k 2 /4' 
due to the Bohm potential, which only becomes important for short wavelengths, 
instability condition then simply becomes the ordinary Jeans instability condition, 



m„TO, m 



For parallel propagation, the 



(17) 



Formally, this means that our system is always unstable for sufficiently long wavelengths. However, we may interpret 
the condition (fl7|) such as to give a required size of the system for an instability of this type to be possible. In 
particular, for the above calculation to be applicable, the system must be approximately homogeneous over length 
scales of the order of the Jeans length, Aj = 27r/fcj, where kj = (AnGpo/c 2 ) 1 ' 2 . Next we consider orthogonal 
propagation, the instability condition becomes, 



V 



1 



— C\ < 0. 



i + x 

In this case the system may be unstable even if the gravitational effect is omitted, provided the condition 

„2 



1 



+ X 2 



1 



C\ 



(18) 



(19) 



is fulfilled. The condition p9[) then suggests that the plasma is unstable for all length scales. In particular even in the 
short wavelength regime, where our calculation is of special interest if we assume that the plasma system of study has 
a limited extension. For short wavelengths, on the other hand, we expect the Bohm potential to be of importance. 
Noting from (|16|) and (fl8|) that the Bohm potential acts stabilizing for sufficiently short wavelengths, this term is 
therefore kept below. 

Next we normalize the parameters according to 



kjc s 



k 



Ca = Ca/c s 



Q 



n 2 k 2 



4m e miC 2 ' 



4 



such that the dispersion relation can be written as 



-4 -2 
UJ — UJ 
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The growth rate 7 = Im Q as a function of k z and has been illustrated in Fig. 1 for (a) the classical Jeans 
instability and (b) the Jeans instability modified by the presence of an external magnetic field. The external magnetic 
field will decrease the instability in the direction orthogonal to the magnetic field. This will generate an oblate 
spheroidal collapse. This is expected since the magnetic forces are anisotropic so the magnetic pressure primarily 
work in directions perpendicular to B. While the oblate shape of the plasma cloud is an effect that occur in the later 
(nonlinear) stage of the collapse, we note that the stabilizing influence of the magnetic field can be seen already in the 
initial linear stage of the instability. In particular the dependence of the growth rate on the wavenumber shows that 
perturbations perpendicular to the magnetic field is less unstable, as is seen by comparing Fig. (la) and Fig. (lb). 
Furthermore, we note that the electron non-locality effects described by he Bohm potential has an over all stabilizing 
effect of the system which is illustrated in Fig. (lc). 

Effects due to the electron spin has a destabilizing effect on the system, in particular in directions almost, but not 
quite, perpendicular to the magnetic field. This instability stems from the magnetic attraction of individual spins (i.e. 
magnetic dipole moments) and has been investigated in Ref. [3f| . The effect of electron spin on the Jeans instability is 
demonstrated in Fig. (Id). In this parameter regime the spin magnetization is significant, but the plasma would still 
be stable if gravitational effects were disregarded. This is in contrast to Fig. (le) in which the spin contribution alone 
is sufficient for collapse. Here, stability is regained for short wavelengths through the influence of the Bohm potential. 
For comparison, the spin instability alone, i.e. without the gravitational Jeans contribution, has been illustrated in 
Fig. (1/) with the same parameter values as in Fig. (le). 

The maximum value of |k| = /c max for which we have instability occur at an angle to k^, as can be seen in e.g. 
Fig. (Id). Our investigation is of particular interest when k max > kj. Especially when considering a system of finite 
size, with inhomogencity scale length shorter than the Jeans length, i.e. for Xj (Vno/no) larger than unity. In this 
regime there will typically be no collapse due to the Jeans instability alone, but as shown from the above analysis a 
collapse for short wavelengths (which is well described by the homogeneous background model) can still occur when 
gravitational and magnetization effects are combined. 



In the present study the effect of magnetization due to the electron spin has been combined with Newtonian 
gravity, in order to describe the stability properties of an initially homogeneous plasma. We note that a more 
complete treatment of the Jeans instability would require taking cosmological effects into account j4lj . As expected, 
the gravitational effects tend to dominate on the large scales, i.e. for long wavelengths. For systems of a limited 
size, background inhomogeneities not included in our model can stabilize the longest wavelengths, in which case 
more significance is given to the short wavelength properties of the dispersion relation. Furthermore, for a system of 
moderate temperature and high density, it is seen that the magnetization of the plasma can contribute significantly 
to the instability in the short wavelength regime. For short wavelengths, we also note that the Bohm potential can be 
of importance, providing a stabilizing influence. A case where the present study is of special interest is when the long 
wavelength perturbations are stabilized by inhomogeneities, and the magnetization of the plasma is significant but not 
sufficient to cause an instability by itself. For such a plasma, the combined effects of gravitation and magnetization 
is needed to fulfill the instability condition. 

In the present model we have not included the effect of the Hall current [I^] , or other finite Larmor radius effects [43[ 
that is often introduced to improve the ideal MHD equations. We note that such modifications introduce dispersive 
properties of the waves, similar to the effects of the Bohm potential. However, as far as the stability properties is 
concerned, the moderate temperature high density regime is of most interest, since it is in this regime that the spin 
effects can give a significant contribution to the instability, as seen above. Furthermore, for such parameters the 
effects of the Bohm potential is more important than those of the Hall current. 

Naturally there are several processes that can be expected in a physical collapse that are not captured within the 
low-frequency MHD limit. For instance, electrostatic repulsion, that could have a stabilizing effect on the system, is 
not included due to the assumption of quasineutrality. Furthermore, since our ideal MHD equations have an infinite 
conductivity, magnetic field lines will be frozen-in during collapse. In reality, there will be particle collisions within the 
plasma and thus the conductivity will always remain finite. The finite conductivity will prevent strict flux-freezing, 
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and as a result the magnetic field will not grow as fast as MHD predicts. This will also reduce the effects of spin 
magnetization somewhat. 

As noted above, magnetization effects as well as electron tunneling and exclusion principle effects, are particularly 
important in low temperature high density systems. In particular only electrons with energies larger than the Fermi 
energy will be free to contribute to the magnetization process. Thus, the magnetization will need to contain an 
effective density of electrons, reduced by the number of electron below the Fermi surface. As such, a different effective 
scaling of the magnetization as a function of density will be obtained. Moreover, correlation effects in dense plasmas 
is of importance, as well as relativistic effects. These are examples of possible modifications that could be introduced 
in modeling of such quantum plasmas in the future. 

The analysis here has been performed through a linearization of the perturbed parameters. More sophisticated 
methods exists which can give qualitative descriptions of the non-linear dynamics of a collapse. The use of integral 
virial theorems, for instance, has proven to often give a more reliable global view of the interaction of magnetic and 
gravitational fields in the non- linear domain [44| . Still, for a more complete understanding of a physically realistic 
collapse, the use of advanced computer simulations is often an necessity. 
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(e) (f) 

FIG. 1: The growthrate has been plotted for (a) the classical Jeans instability and (b) the Jeans instability modified by the 
presence of an external magnetic field (C'a = 2 in Fig. (b) — (/)). In (c) we have added non- locality effects of the electrons 
described by the Bohm potential (Q — 0.5), and (d) includes electron spin effects in the plasma description (a = 0.4). In 
(e) both spin and non-locality effects are included in the model in a regime (a = 0.75, Q = 0.5) where the plasma would be 
unstable even without graviational effects, Fig. (/). 



